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Abstract: In (3 + l)-dimensional SU{Nc) Yang-Mills (YM) theory, the Chern- 
Simons diffusion rate, Tcsi is determined by the zero-momentum, zero-frequency 
limit of the retarded two-point function of the CP-odd operator tr [F A F] , with F 
the YM field strength. The Chern-Simons diffusion rate is a crucial ingredient for 
many CP-odd phenomena, including the chiral magnetic effect in the quark-gluon 
plasma. We compute Vcs in the high-temperature, deconfined phase of Improved 
Holographic QCD, a refined holographic model for large-iVc YM theory. Our result 
for Tcs/{sT), where s is entropy density and T is temperature, varies slowly at high 
T and increases monotonically as T approaches the transition temperature from 
above. We also study the retarded two-point function of tr [F A F] with non-zero 
frequency and momentum. Our results suggest that the CP-odd phenomena that 
may potentially occur in heavy ion collisions could be controlled by an excitation 
with energy on the order of the lightest axial glueball mass. 
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1. Introduction 

(3+l)-dimensional SU{Nc) Yang-Mills (YM) theory has an infinite number of de- 
generate classical vacua distinguished by a topological invariant, the Chern-Simons 
number, Ncs- Normalizing the YM kinetic term as — ^tr[F^j,F'^'^], Nqs is 



^cs = j ^^i*^ 



where i,j,k = 1,2,3 and the trace is over gauge indices. A change in the Chern- 
Simons number is thus 

AAcs = [ d^xqix"), (1.2a) 



q{xn ^ ^tr [i^ A F] = ^ e--trF,.F,., (1.2b) 

where = {t,x). In a state invariant under translations in space and time, the 
rate of change of Acs per unit volume V per unit time t is called the Chern-Simons 
diffusion rate, denoted Fes, 



, ^((AAcsT) 
Vt 



d'^x (g(x'^)g(0))w, (1.3) 



where the subscript W denotes the Wightman function. In an equilibrium state with 
non-zero temperature T, let Gr{uj, k) denote the retarded Green's function of q{x'^) 
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in Fourier space, with frequency u and spatial momentum k. In such states, eq. (|1.3| ) 
can be rewritten as 

2T 

Tcs = -hm — lmGR(uj,k = 0). (1.4) 

Gauge field configurations for which J d'^xq{x^) is non-zero produce a non-zero 
AA^'cs- At zero temperature, such gauge field configurations, called instantons, rep- 
resent quantum tunneling events between vacua. At both zero and non-zero T, the 
contribution of instantons to Fes is exponentially suppressed [1,2]. When T is non- 
zero, however, classical thermal fiuctuations can also produce a non-zero AA^^cs, for 
example by exciting unstable gauge field configurations called sphalerons [3,4] which 
generate non-zero AA'cs upon decay. Such classical thermal processes are not ex- 
ponentially suppressed [5-7]: in YM perturbation theory Fes oc log(A4) T'^, where 
\t = g'^Nc is the 't Hooft coupling [8-11]. 

In YM coupled to fundamental-representation fermions, J d^x q{x^) also con- 
tributes to chiral anomalies in global symmetries. In the electroweak theory, gauge 
field configurations with non-zero j d'^x q{x^) play a role in electroweak baryogene- 
sis [12,13], while in Quantum Chromodynamics (QCD), for sufficiently high T they 
may play a role in generating bubbles of net chirality (more left-handed than right- 
handed quarks, for example), in which parity, P, and charge conjugation times parity, 
CP, are broken [14].^ 

Such CP-odd domains in hot QCD may have observable consequences in heavy 
ion collisions at the Relativistic Heavy Ion Collider (RHIC) and Large Hadron Col- 
lider (LHC). These collisions produce a hot soup of QCD matter with T on the order 
of two to four times the QCD crossover temperature. The resulting state behaves as a 
nearly-ideal fiuid of strongly-interacting quarks and gluons, the quark-gluon plasma 
(QGP) [15-17]. A non-central collision may produce a QGP with non-zero angular 
momentum and hence a magnetic field, both pointing perpendicular to the reaction 
plane (spanned by the beam axis and impact parameter). In the presence of a mag- 
netic field, a net chirality will produce an electric current parallel to the magnetic 
field, due to the axial anomaly. This is the Chiral Magnetic Effect (CME) [18, 19]. 
A detection of the CME in heavy ion collisions would thus be a detection of CP-odd 
processes in QCD. 

One observable consequence of the CME in a heavy ion collision is charge sepa- 
ration: positive charges will move to one side of the reaction plane, negative charges 
to the other. We know from experiment that the strong interactions preserve P and 
CP, however, so any charge separation from CP-odd sources will, over many events, 
average to zero. An observable that could serve as a "smoking gun" for the CME is 
thus hard to find. For heavy ion collsions at RHIC and LHC the focus so far has 
been on three-particle correlations [20-22], which indeed indicate that charge separa- 

^ Sometimes Fes is also called the "sphaleron transition rate" or, in the context of electroweak 
baryogenesis, the "baryon number violation rate." 
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tion occurs in heavy ion collisions. These correlations are sensitive to event-by-event 
charge separation from both CP-odd and CP-even processes, however, making a 
positive identification of a signal from the CME difficult [23]. In short, to date the 
experimental evidence for the detection of the CME in heavy ion collisions at RHIC 
and LHC is inconclusive. 

The experimental situation raises a number of urgent questions for theorists. 
Can we compute the size of the signal from the CME, relative to backgrounds? How 
will that signal depend on temperature, magnetic field, centrality, etc.? Clearly an 
auxiliary question is: how big is the rate of chirality production, which is oc Fes, in 
a heavy ion collision? 

Unfortunately, Fes is difficult to calculate for the QGP, for the same reasons that 
the shear viscosity, rj, is difficult to calculate. The quarks and gluons are strongly- 
interacting, so perturbation theory is a priori unreliable. Calculating transport coef- 
ficients, such as Fes and rj, from lattice QCD requires a problematic analytic contin- 
uation from Euclidean signature.^ Currently no reliable technique exists to compute 
Fes or V for QCD at the temperatures reached in the QGP. 

An alternative approach is holography [26-28], which equates certain strongly- 
coupled gauge theories in the large- A'^c limit with weakly-coupled theories of gravity in 
spacetimes of one higher spatial dimension. A deconfined thermal state of the gauge 
theory is dual to a black hole spacetime [29] , and transport coefficients are relatively 
straightforward to calculate [30-32]. Remarkably, the ratio of ri to entropy density, 
s, for any theory dual to higher-dimensional Einstein gravity is rj/s = l/(47r) [33], 
which is close to the estimate for rj/s for the QGP extracted from data [34]. Such 
universality serves as encouragement for computing other transport coefficients, like 
Fes, from holography. 

Previous calculations of Fes in holography employed "top-down" models, i.e. 
models descending from a known string theory or supergravity construction. The 
best-understood example is A/" = 4 supersymmetric YM (SYM) with large Nc and 
large Xt, dual to supergravity in an Anti-de Sitter (AdS) space, for which Fes oc 
T"^ [30] . Other holographic calculations included the effects on Fes due to a mag- 
netic field [35] or confinement [36]. To our knowledge, in all previous cases the 
holographic results for Fes were ultimately fixed by some underlying (perhaps "hid- 
den" [36]) conformal symmetry. 

In this paper we compute Fes in Improved Holographic QCD (IHQCD) [37-44], 
a holographic model of large- A'c YM theory. The model is "bottom-up," i.e. does not 
descend from a known string theory or supergravity construction, but is tailored to 
model string theory systems very closely, unlike other bottom-up models. The bulk 

^In fact, Fes may be more difficult to calculate from lattice QCD than other transport coeffi- 
cients. The operator q{x'^) obeys various constraints. For example, ANqs ~ J d'^xq{x'^) must be 
an integer, the second Chern character. Defining a lattice version of the operator q{x^) that obeys 
all of the constraints can be difficult, as discussed for example in refs. [11,24,25]. 
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theory is Einstein-dilaton gravity, where the dilaton $ is dual to trF^j^F^^. A non- 
trivial dilaton solution will describe non-trivial running of the YM coupling, hence 
the choice of dilaton potential is crucial. The simplest choice involves only two free 
parameters, which can be adjusted such that the model reproduces both the T = 
glueball spectrum and the thermodynamics of large-Nc YM, including a first-order 
deconfinement transition at a critical temperature T^. In particular, the model has 
no (hidden) conformal symmetry. We briefly review IHQCD in section ^. 

In IHQCD the operator q{x'^) is dual holographically to an axion field in the 
bulk [37,38,40,42,44]. Defining for convenience a holographic 't Hooft coupling 
A = e*, the normalization of the axion's kinetic term includes a dilaton-dependent 
factor, Z{X). In principle, Z{X) could be fixed by matching to lattice results for the 
Euclidean correlator of q{x^), as we explain in section 0. We work instead with several 
simple choices for Z{X), in part to study the generic behavior of Fes in holographic 
models. Specifically, we consider a Z{X) with two free parameters, which we fix 
by demanding that the model match large- A^^^ YM lattice results for the topological 
susceptibility and for axial glueball mass ratios to within one sigma. 

In section || we compute Fes in the high-temperature, deconfined phase of 
IHQCD. Letting s denote the entropy density and Xh the value of A at the black 
hole horizon, our result for Fes is of the form 

rc. = ^f ^(A.). (1.5) 

Figs. H, D, and |] show our numerical results for Tcs/{sT/N^) = Z{Xh)/{2-n). For 
our choices of Z{X), the value of Z{Xh) is bounded from below as a function of T by 
its value in the T — )■ oo limit, and increases monotonically as T approaches Tc from 
above, with most of the increase occurring between 2Tc and T^. We will argue that 
such behavior is generic in a large class of confining theories with classical gravity 
duals. In a scan through various choices of Z{X), each of which reproduces the first 
two axial glueball mass ratios to within one sigma, we find that the increase can be 
as large as 60%. For our optimal choice of Z{X), which provides the best fit to the 
lattice results for the first two axial glueball mass ratios, the increase is only 0.01%. 

To obtain Fes, we compute the low-frequency limit of Gji{uj, k = 0) holograph- 
ically. In section ^ we initiate the study of GR{u,k) at non-zero u and \k\, in the 
T > Tc regime. We focus in particular on lmGii{uj,k), which is proportional to 
the spectral function of q{x^). After suitably subtracting the high-frequency asymp- 
totics, by computing the difference in the value of the correlator at two temperatures, 
our results suggest the presence of a reasonably long-lived excitation with energy on 
order of the lightest axial glueball mass at T = 0. That is sufficiently light to prompt 
the speculation that perhaps such an excitation could dominate many CP-odd phe- 
nomena in the QGP created in heavy ion collisions. 

In section ^ we summarize our results and discuss directions for future research. 
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2. Improved Holographic QCD 



The holographic model that we consider as the dual to pure large- A'^c YM is (4+1)- 
dimensional Einstein-dilaton gravity with a well-chosen dilaton potential [37-44]. In 
terms of the holographic 't Hooft coupling A = e*, the bulk action is 



S = MlNl J d^x^/^ 



2 

+ 'S'bdry, (2.1) 



where Mp is the Planck Mass, related to the (4+l)-dimensional Newton's constant 
Gs as Mp = l/(167rG'5A^^), g and R are the determinant and Ricci scalar of the 
bulk metric, ^(A) is the dilaton potential, and S'bdry represents all boundary terms, 
including the Gibbons-Hawking term as well as counterterms. 

If V^(A) = 12/£^ with a constant length scale then the equations of motion 
arising from eq. ( |2.1| ) admit a solution with constant A and an AdS metric with 
radius of curvature £, 

ds\^s = ^ (^^^ - d'>^^ + d^"^) > < r < cx). (2.2) 

Here r is the holographic radial coordinate, dual to the field theory energy scale: the 
region near the AdS boundary at r — )• is dual to the ultra-violet (UV) of the field 
theory, while the region near the Poincare horizon at r — >■ oo is dual to the infra-red 
(IR). Such a solution describes a conformal field theory. 

For non-trivial V^(A), the equations of motion admit vacuum solutions in which 
A depends only on r and the metric takes the form 

ds'^ = bo{rf{dr^ - dt^ + dx^), < r < oo, (2.3) 

with warp factor hQ{r). In IHQCD we demand that as r — )■ the metric approach 
that of AdS, hQ{r) — )■ ^/r (up to corrections logarithmic in r) and that A vanish 
logarithmically, A — )■ — 1/ logr, to mimic the running of the large- A'^c YM coupling. 

Large- Ai'c YM approaches a free theory in the UV, so we expect the holographic 
dual in the r — )■ region to be a string theory, not just a classical gravity theory 
like IHQCD. On the other hand, in large- A^^c YM, A^ diverges in the IR, so a classical 
gravity theory may be a reliable description in the r — )■ oo region. IHQCD is intended 
to be such a low-energy effective description of large- A^^c YM, reliable in the r — )■ oo 
region. In practice, the role of the r — )■ region in IHQCD is simply to provide 
boundary conditions for the fields in the r — >■ oo region. We impose those boundary 
conditions at a cutoff, i.e. at some small but finite r = e. We then compute low-energy 
quantities that are insensitive to the cutoff, some of which we match to large- A"c YM, 
while the rest are predictions of the model. A more detailed discussion of these issues 
appears for example in ref. [41]. 
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Using classical gravity in the r — )■ region has an important consequence, how- 
ever: IHQCD will actually be dual to a theory that flows to a non-trivial UV fixed 
point. Generically, the UV physics of large- A^c YM and IHQCD will thus be different. 
For example, in IHQCD, rj/s = l/(47r) [33], which is much smaller than the high-T 
perturbative result for rj/s in large- iVc YM [45]. Nevertheless, in order to match 
IHQCD to known results for IR quantities in large-A'c YM, we must match some 
quantities in the UV. For example, to reproduce lattice results for the free energy 
of large- iVc YM for T > Tc with the correct normalization, we must demand that 
at high T the free energy of IHQCD obey a Stefan-Boltzmann law. That require- 
ment fixes the value of i in the asymptotic AdS region in units of the Planck mass: 
(M/)-3 = 457r2 [39,40]. 

By matching to another UV quantity, the perturbative large- A^^c YM /3-function, 
we can also constrain V^(A). In the r ^ region, where A is small, ^(A) has a regular 
series expansion 

12 

V{X) = -{l+voX + viX^ + O (A^)) . (2.4) 

Committing to an identification of the field theory renormalization scale E = Eq bQ{r), 
where Eq can be fixed by matching to the lowest glueball mass or to the result for 
Tc from lattice large-A^c YM, we can fix the coefficients Vq and Vi in terms of the 
coefficients of the perturbative large- A^c YM /3-function [37,38,44]: 

22 ^1 

/3(A*) = -/3oA?-/3iA? + 0(A^), /^o = '^^ = m^O' (2.5a) 

^0 = IPo, V, = + (2.5b) 

In the vacuum solutions, generically A diverges as r — )■ oo. The large-A expansion 
of V^(A) must take the form V^(A) oc As i/log A in order for the glueball spectrum to 
be gapped and discrete with asymptotically linear trajectories [37,38,41,44]. With 
this asymptotic form for V{X), as r — )■ oo the warp factor and A(r) take the form 

6o(r) oce-^"^/^)', A(r) oc y et^^/^)', (r ^ oo) (2.6) 

where L is a length scale determined by the value of A at r = e. The form of bo{r) in 



eq. ( p^) is sufficient to guarantee that the dual field theory is confining [37,38]. The 
metric actually has a mild singularity^ at r = oo that can be cloaked by a regular 
horizon and hence is a "good" singularity [46]. Moreover, the singularity is repul- 
sive [38,40], which guarantees that the low-energy spectrum and other observables 
are insensitive to the details of the resolution of the singularity. 



■^On the other hand, m the string frame, where the metric scale factor is A^/"^(r)6o(r'), the 
curvature approaches zero as r — >■ cx). 
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The black hole solutions of the model defined by the action in eq. (|2.1| ) have 
non-trivial A(r) and a metric of the form [39,40] 



ds^ = b{ry 



fir) 



f{r)dt' + dx' 



< r < Th. 



(2.7) 



The surface r = r/j is the horizon, where f{rh) = 0, and the corresponding Hawking 
temperature is T = 47r/'(r/i). Black hole solutions only exist for temperatures above 
a value Tmin, and in fact two branches of solutions exist, the large and small black 
holes (comparing r/j to i). Fig. |] depicts the typical form of T as a function of 
Vh, including the two branches of black hole solutions. For both large and small 
black holes, as r — )■ r^, the warp factor b{r) asymptotes to a constant whose value 
determines the entropy density, s = b{rhY / (AG^), and as r — )■ 0, 6(r) — > r/i, up to 
0{r^) (times logarithmic) corrections, indicating that in the field theory the thermal 
energy density and pressure are both of order N^. 

T 



min 



IBig black holes 


Small black Holes 


L 









min 



Figure 1: Schematic plot for the typical form of the black hole Hawking temperature T as 
a function of the horizon position r^, for a generic choice of V{X) (with the correct small- 
and large- A asymptotics) . The temperature exhibits a minimum, T^nm? '^t t^i^i^ which 
separates the large black hole (r/j < rmin) from the small black hole (r/j > rmin) branches. 

In large black hole solutions, A(r) decreases monotonically as T increases, so that 
A — )■ as T — )■ oo. In small black hole solutions, A(r) increases as T increases. In 
particular, as discussed in refs. [40,42], the value of A(r) at the horizon, = A(r/j), 
is a monotonically increasing function of r^, so a plot of T versus Xh is qualitatively 
similar to fig. ||: in the T — )■ oo limit, A^ — ?■ on the large black hole branch and 
A/i — 7- oo on the small black hole branch (see for example fig. 2 (a) of ref. [40]). 
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If we Wick-rotate to a compact Euclidean time direction of length 1/T, then for 
T > Tmin three bulk solutions exist: the Wick-rotated version of eq. (|2.3| ), which 
describes a thermal gas of gravitons and is dual to a confined state, and the Wick- 
rotated large and small black holes, which are dual to deconfined states. To determine 
which solution is thermodynamically preferred at any given T, we must determine 
which has the smallest on-shell Euclidean action, dual to the field theory's free en- 
ergy (times 1/T). As shown in refs. [39,40,42], the small black hole solutions are 
never thermodynamically preferred, but at some Tc > Tmin the large black hole solu- 
tions become thermodynamically preferred. Indeed, the system exhibits a first-order 
Hawking-Page type transition at Tc, dual to a confinement-deconfinement transition. 

In general, for a given potential V^(A) we cannot solve the equations of motion 
arising from eq. (|2.1| ) exactly, so we resort to numerics. Here we will only sketch 
our numerical procedure, which is described in detail for example in ref. [42]. At 
the cutoff r = e we impose a Dirichlet condition on each field, and in particular we 
demand that the metric take the AdS form. We then fix the remaining integration 
constants, including A^, by a shooting algorithm. Given a choice of V{X) and the 
Dirichlet conditions at r = e, we obtain a one-parameter family of solutions labeled 
by T, or equivalently by Xh- Following refs. [42,43], in our numerics we use a simple 
form for V^(A) with the correct small- and large-A asymptotics. 



I + V0X + ViA^/yiog (1 + V2A4/3 + 1/3 A2) 



(2i 



Expanding eq. ( pTq ) about A = and matching to eq. ([2^^), we find vq = Vq and 
f 1 = Via/T^- The coefficients Vq and V2 can be determined in terms of Vi by imposing 
the conditions in eq. ( |2.5b| ). The potential thus has two free parameters, Vi and V3. 



We fit these two parameters by matching to lattice results for two thermodynamic 
quantities in large- iVc YM: the latent heat of the deconfinement transition, which 
is proportional to the entropy density at the transition, s(Tc) / (N^Tj^) ~ 0.31 [47], 
and the pressure at T = 2Tc [47-49]. Upon fixing Vi and V3 in this fashion, IHQCD 
describes very well both the T = glueball spectra (0+"*" and 2++) as well as the 
finite T thermodynamics of laige-Nc YM [42,50]. 

The instanton number density operator, q{x^) in eq. (|1.2b|) , is dual to a bulk 



pseudoscalar, the axion a (as in many top-down models). In the field theory, the 
source for q{x'^) is an angular variable, the ^-angle. As a result, the action of the 
bulk axion must be invariant under shifts of a, and hence must depend only on 
derivatives^ da. General arguments in string theory and in YM theory, including the 
argument that 6 dependence should appear in the YM vacuum energy only at order 
one in the large- A^c limit rather than at order A^^ [51], imply that the axion action 



^Instanton effects may produce a non-trivial axion potential, such as a term cos a. These in- 
stanton effects are exponentially suppressed in the large- iVc limit, however. 
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is suppressed by Oil/N"^) compared to the action S in eq. (pTTI ) [37,38,41,44]. We 
thus add to the model an axion with an action Sa of the form [37,38,41,44] 



(2.9) 



where, following the rules of effective field theory, we have included a dimensionless, 
A-dependent normalization function, Z{X), consistent with the symmetries. 

Being a massless pseudo-scalar, in an expansion of a(r) about r = 0, the leading, 
non-normalizable term is a constant, which is proportional to the YM ^-angle defined 
in the UV, 



where in top-down models the proportionality constant k will be fixed, but not in 
bottom-up models. In other words, in our model the normalization of the operator 
dual to a is ambiguous: a is dual to q{x^)/ n. Nevertheless, by fixing the normaliza- 
tion of the topological susceptibility we will be able to compute two-point functions 
of q{x^) unambiguously, as we explain below. 

To specify Sa completely we must specify Z{X). In principle, Z{\) can be fixed 
as follows. First, perform a lattice calculation of the Euclidean two-point function 
of g(x^) with non-zero T for some set of frequencies. Second, compute the same 
Euclidean two-point function holographically for all frequencies for some choice of 
Z{\). A least squares fit of the holographic results to the lattice results should then 
determine Z{X). To study generic behavior of holographic models, we will instead 
proceed by using simple forms for Z{X) that we constrain by matching to lattice 
results for the topological susceptibility and axial glueball mass spectrum. Notice 
that matching to any lattice data will always have room for improvement: lattice 
definitions of q{x^) generically suffer from power-law divergences that dominate in 
the continuum limit, making lattice calculations of correlators of q{x^) noisy [24]. 
Accurate calculations may be possible in the near future^. 

We can constrain Z{X) as follows. Since Z{X) is the coefficient of a kinetic 
term, we demand that Z{\) > 0. We can also constrain Z(A)'s small- and large-A 
asymptotics [37,38,41,44]: 



where Zq is a dimensionless constant. The small- A form follows from the rules of 
effective field theory: a constant is the most general allowed term. The large A 
behavior is fixed by glueball universality [38] . Various towers of glueballs have linear 
asymptotic trajectories: for large excitation number n, their squared masses go as 

^We thank F. Bruckmann, H. Panagopoulos, and A. Schafer for discussions on this issue. 



a(r = 0) 



(2.10) 




A — > oo. 



(2.11) 
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(m^)^ = c*n + ■ ■ ■ , with constants c*, where the integer i labels different towers. 
Glueball universality is the statement that all the slopes d are similar, i.e. do not 
depend on i. That is automatic for the O"*""*" and 2~^~^ glueballs. Requiring the same 
for the O""*" glueballs forces Z{X) to go as A"^ at large A [38]. 

We will use the simplest form of Z{X), also used for example in ref. [42], 

Z(A) = Zo(l + C4A^), (2.12) 

where C4 is a dimensionless constant. To fix Zq we match to the large- iVc YM lattice 
result for the Euclidean topological susceptibility, defined in terms of the T = 
vacuum energy density £{9) as 

X^^-§^ = J d'x{q{xnqm^. (2.13) 

where the subscript E denotes the Euclidean correlator. The holographic result for 
X is'' [38], 

Jo 6f,(r)Z{A{r)) 

Clearly x will be proportional to k^Zq. Thus, for any given value of the parameter C4, 
matching the holographic result for x to the lattice result, x ~ (191 Me V)"^ [24,52], 
fixes the product k'^Zq. On the other hand, since the locations of poles in the two- 
point function of g(x'^) are independent of the overall normalization k'^Zq, we can 
fix C4 independently by matching the mass of the lowest 0""*" glueball to the lattice 
result of ref. [53], 

mo-+/mo++ = 1.50(4). (2.15) 

The resulting values are^ [42] 

k'^Zo = 33.25, C4 = 0.26. (2.16) 

These values can then be used to predict the masses in the full tower of 0""*" glueballs. 
As shown in refs. [37,38,41,44], the holographic result for the first excited 0""*" glueball 
mass, mo*-+, agrees very well with the lattice result [53], 

mo.-+/mo++ = 2.11(6). (2.17) 

Crucially, notice that by fixing the normalization of x we have fixed the normal- 
ization of any two-point function of q{x'^), and thus have eliminated the normalization 
ambiguity mentioned below eq. ( |2.10| ). In other words, the holographic calculation 



^The holographic calculation of x in ref. [38] assumed k = 1. Here we allow for arbitrary k. 
''The result for k^Zq in ref. [42] was too large by a factor of four, producing an erroneous result, 



K Zq = 133. In eq. (2.16) we present the correct value, k Zq ~ 133/4 = 33. 2£ 
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of the two-point functions of q{x^) will only depend on the combination k'^Zq (as we 
will see explcitly in section |^), which we have fixed to the value in eq. (|2.16|) . 

Solutions for a as a function of T were studied in refs. [40,42]. When T = 0, 
a non-trivial UV ^-angle forces a{r) to be non-trivial. The resulting normalizable 
solution then indicates that the non-zero UV 6'-angle flows to zero in the IR, as 
shown in fig. H, and additionally triggers a non-zero The T = solution for 

a{r) is unchanged when T < T^: Wick-rotating the metric in eq. ( p.3|) to a compact 
Euclidean time does not affect the static solution a{r). Such behavior is expected 
in a confined phase at leading order in Nc, due to large-A^^c volume independence. 
When T > Tc, however, the only non-singular solution for the axion is a constant, 
a{r) = kO, indicating that {q{x)) = 0, in agreement with evidence from lattice data 
for large- A/'c YM [24]. 




100 200 300 400 500 600 



E (MeV) 



Figure 2: The normalizable solution a(r) for the axion at T = 0, expressed as a running 
^-angle normalized to the UV value, as a function of the energy scale, E{r) = EQb{r). We 
fix Eq by matching our holographic result for Tc to the large- A'^c YM lattice result. 

What is the behavior of the topological susceptibility as a function of tempera- 
ture, x(T)7 When T < Tc, xiT) is independent of T, i.e. takes the same value as at 
T = 0, eq. ( p.l4| ), again due to large-Yc volume independence. When T > Tc, the 



holographic result for the topological susceptibility is 



xiT) = .r, {T > Tc) (2.18) 

Jo b^r)f{r)Z{X{r)) 

The denominator on the right-hand-side of eq. ( p.l8| ) diverges at the black hole 
horizon, so in fact x{T) = when T > Tc, up to 0(e~^=) corrections [40,42]. 

We will also consider a form for Z{X) more general than that of eq. ( p.l2| ). On 
the large black hole branch, if T is large then A is small, in which case we expect 
the largest poljTiomial correction to the Z{X) in eq. ( |2.12| ) to be a term linear in A, 
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hence we consider 

Z(A) = Zo (1 + ciA + C4A^) , (2.19) 

where ci is a dimensionless constant, which we choose to be positive. If we continue 
to fit only to the lattice result for the lowest glueball mass, we find a substantial 
degeneracy (which is not surprising, given that we have introduced an additional 
parameter, ci). Specifically, for any positive value of Ci, a value of C4 exists such 
that, upon matching to the lowest axial glueball mass in eq. (|2.15|), the value of the 



first excited axial glueball mass is in rough agreement with the value in eq. ( |2.17| ), 
exhibiting at most a 3% discrepancy, as shown in fig. ^. To constrain ci we will 
demand that our holographic results for the axial glueball masses fall within one 
sigma of the lattice values in eqs. ( p.l5|) and (|2.17|) . That results in the constraints 



< ci < 5, 0.06 < C4 < 50. (2.20) 

In fact, the optimal values, which provide the best fit, are the ones in eqs. ( p.l2| ) 
and {^A§: (01,04) = (0,0.26). 

As we have seen, the function Z{X) must be non-negative and is constrained in the 
A — and A — )■ 00 limits. For intermediate values of A, the most natural assumption 
is that Z{X) is monotonic. At least, we are not aware of any compelling evidence 
for the existence of maxima or minima in Z{X). Our choices for Z{X) were thus 
monotonic functions of A, namely polynomials in A with strictly positive coefficients. 
To test the effect of maxima and minima in Z{X), we considered two changes to the 
Z{X) in eq. (|2.19|) . First, we allowed slightly negative ci, while maintaining Z{X) > 0. 



Second, we introduced a maximum by a adding a Gaussian peak to Z{X). In each 
case we computed the axial glueball mass spectrum. After matching to the lattice 
result for the lowest axial glueball mass, we found that the fit to the first excited 
axial glueball mass was worse, deviating from the lattice result by about 10%. We 
consider that a preliminary indication that monotonic Z{X) may indeed be the best 
choice. We leave more thorough tests for future research. 

Our assumption that Z{X) is monotonic in A determines the qualitative behavior 
of Z{X) as a function of T. On the large black hole branch, as T — t- 00, A — t- 0, and 
as T decreases towards Tc, X increases monotonically. As a result, for our choices of 
Z{X) — simple polynomials in A with positive coefficients — when T — > 00, Z{X) — t- Zq, 
and when T — )■ Z{X) will increase monotonically. As functions of T, our Z{X) 
are thus bounded from below by their value in the T — )■ 00 limit: Z{X) > Zq. The 
behavior of Z{X) as a function of T will translate directly into the behavior of Fes as 
a function of T, as we will show in the next section. In particular, the dimensionless 
combination Fcs/('5T) will be bounded from below by its value in the T — )■ 00 limit, 
and will increase as T — T^. In the next section we will also present a more general 
argument that Fcs/('5T) must increase as T approaches from above. 
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Figure 3: Our holographic results for the masses of the glueball states with 
excitation number n, normalized to the lowest 0"*"*" glueball mass, obtained by vary- 
ing the coefficients ci and C4 in the Z{X) in eq. ( ^.19 ). From the top (red) 



dots (visible only for n = 2 and n = 4) to the bottom (blue) dots, (01,04) = 
(0,0.26), (0.5, 0.87), (1,2.2), (5, 24), (10, 75), (20, 230), (40, 600). The lowest axial gluebah 
mass, n = 1, is always fixed to be the value in eq. ( |2.15| ). The two horizontal blue lines 
with surrounding blue bands indicate the results and errors, respectively, of the large-A^'c 
YM lattice calculations for the masses of the lowest and first excited states, n = 1 and 



n = 2 (see eqs. (2.15) and (2.1?)) [53]. Only the mass of the n = 2 state is appreciably 



sensitive to changes of ci and C4, differing from the lattice result by 3% at most. 

3. The C hern- Simons Diffusion Rate 



We will compute Tqs using eq. (|1.4|), rewritten as 



2T 

Tcs = -K^ lim — lmGR(u,k = 0), (3.1) 

where GR{uj,k) is the retarded two-point function of q{x^)/K, the operator dual to 
our axion a. 

In holography, the on-shell bulk action is the generating functional for field theory 
correlation functions [27,28]. To compute the two-point function Gr^lo,}:) in the 
high-temperature, deconfined phase of IHQCD, we must solve the linearized equation 
of motion of the axion in the black hole spacetime with metric in eq. (|2.7| ), with 
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T > Tc. We thus introduce a fluctuation of tfie axion, Sa{r,x^), wliere = {t,x). 
Wlien T > Tc, the background solution for the axion is trivial, hence the linearized 
equation of motion for Sa{r,x^) is simply 

dr [Z(A(r))v^/"a,fc(r,x^)] + g^'^d^.d^ 6a{r,x^') = 0, (3.2) 



Z(A(r))v/^ 

where the metric is that of eq. (|2.7| ). Notice in particular that Sa will not couple 
to the fluctuations of any other fields because the background solution preserves CP 
and the axion is the only CP-odd field in the bulk. We must solve eq. ( p.2| ) with 
Dirichlet boundary condition at the asymptotically AdS boundary and with in-going 
wave boundary condition at the horizon [30]. The solution takes the form 

6a{r, xn = e^'^ 5a(r, F) a(F), (3.3) 

— * 

where = (w, k) and where a{k^) is fixed by the Dirichlet boundary condition, 

hm5a{r,xn = J (0ie^'^a(P), (3.4) 

while 6a{r, k^) obeys the equation 
1 



dr [Z{\{r))^g''dM{r,k^)\ - g^'k.k, Sa{r, k^ = 0, (3.5) 



with unit normalization at the asymptotically AdS boundary, limr-^^o 5 a {r,k^) = 1, 
and in-going wave boundary condition at the horizon. The on-shell axion action is 
then 

/ 74 7 
^a(-P)^(r,P)a(P)^ , (3.6) 

where 

J^(r,fc'^) = — ^ 5a{r, -k^") Z{\{r)) g'"'- dr5a{r, k"). (3.7) 



2 

The retarded Green's function is then [30] 

Gji(uj, k) = -2 lim J^(r, P). (3.8) 

r-i>0 



To compute Fes, we need to solve eq. ( |3.5| ) with k = and with small u. We 
will do so in two ways, first using near-horizon matching and second using the mem- 
brane paradigm, following ref. [32]. In each case we can determine Pes analytically, 
essentially because 6a is a massless fluctuation. 

In the near-horizon matching technique, we first solve eq. (|3.5| ) with u = and 
then expand the solution near the horizon. We then reverse the order of operations, 
solving the equation in the near-horizon region and then expanding the solution in 
u. Finally, we match the two solutions to obtain J^{r, k^). 
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When k = and u = the solution of eq. (|3.5|) is 

r dr' 
Sa^C, + C,J^ Z(A(rO)MrOV(r')' 

with constant coefficients Ci and C2. The second term on the right-hand side of 
eq. ( p.9|) diverges as r — )■ r/^. As a result, when w = a normalizable solution must 
have C2 = 0. When u is small but non-zero, a normalizable solution may have 
C2 oc u. Plugging eq. (|3.9| ) into eq. ( p.7|) , we find 

limJ^(r,A;^) = ^CiCs. < T, = 0) (3.10) 

r—^O 2 

We will choose Ci = 1 so that our 6a has unit normalization at the asymptotically 
AdS boundary. Our task is thus to determine C2. Expanding the solution in eq. (|3.9| ) 
around the horizon, we find 

where f'{rh) = AttT. Now we reverse the order of operations. Expanding eq. ( |3.5| ) in 
(r/i — r), we find the solution in the near-horizon region, 

Sa = C+{rh- r)^ + C^{rh - r)~^ , (3.12) 

with coefficients C± that depend on u but not on r. We set C+ = so that the near- 
horizon solution is an in-going wave [30]. Now we expand the solution in eq. ( p.l2| ) 
for small u: 

5a = C^- C_ log(r, - r) + 0{uj^/T^). (3.13) 
By matching the constant and logarithmic terms in eqs. ( |3.11|) and ( p.l3|) , we find 

Ci = C_, C2 = -iujZ{\H)h{rhfC-. (3.14) 

Setting Ci = 1, we obtain limr^Qj^{r,k^^) via eq. ( |3.10D and then Gji{u,k) via 
eq. (IP), 

Gniuj, k = 0) = -ico Ml Z{\r,) b{r^f. {u < T) (3.15) 
We thus obtain our main result for Fes, 



^cs = -k' lim — ImGniu, k = 0) = —^ K^Z{Xh), (3.16) 

w-i-O UJ iV;r 2,71 



where we have used = 1/ (IGnG^N^) and where s = ^jp^ is the entropy density. 



4G5 

.2 ' 



Notice that the normalization of this result is fixed by the product n Zq, which we 
fixed in section by matching to the topological susceptibility at T = 0. 
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The second equivalent, but more efficient, method that we will use to obtain 
Gfi{u,k) is the membrane paradigm [32]. Kubo's formula for the retarded Green's 
function is 

n(w, k) = GR{io, k)6a{io, k), (3.17) 

where Il{u,k) is the one-point function of q{x^)/K in Fourier space. Following 
ref. [32], we extend eq. ( |3.17| ) into the bulk by defining an r-dependent response 
function, 

n(r, k) 



Cir,uj,k) = 



u 6a{r, co, k) 



(3.18) 



where n(r, cu, fc) is the canonical momentum of 5a{r,u,k) with respect to the r- 
foliation of the bulk space-time. 



n(r, cj, k) 



6Sn 



-Ml Z(A(r)) v/=^/^ dM{r^ uj, k). 



(3.19) 



5dr5a 

The retarded Green's function is then proportional to the boundary value of C,: 

GR{uj,k) = -Mluj\imCir,uj,k). (3.20) 



An equation of motion for ( is straightforward to derive using eq. (|3.19|) and 6a' s 
equation of motion, eq. (|3.5|), 



drC 



Z{X{r))^g-r 

UJ 

Z{\{r))h{rYf{r) 



gtt ^Jl 



(3.21) 



To obtain the retarded Green's function Gr{u, k), we must impose regularity at the 
horizon, meaning drC is finite there [32], hence the term in brackets in eq. (|3.21|) must 
vanish^ at r = rt'- 

C(r,) = +zZ(A,)K^h)'. (3.22) 

We can now easily derive Fes- In eq. (|3.21|) we take k = and observe that if a; — ?■ 
then ( becomes independent of r. The value of ( for all r is then the same as the 
value at the horizon, eq. ( |3.22| ), and via eq. ( p.2CI| ) we trivially obtain Gr^oj, k = 0), 
which is identical to eq. ( p.l5| ). We thus find again 



Tcs = lim — ImGR(c<;, A; = 0) = -^^/«^Z(A?,). 



(3.23) 



When k but uj — 0, the boundary condition is modified from that in eq. ( 3.22 ), as discussed 
in ref. [32]. In what fohows, whenever we consider fc ^ we will work with w ^ 0, hence we will 
use the boimdary condition in eq. ( 3.22| ). 
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Our result suggests a natural dimensionless quantity to study, 



sT/N^ 27r ' ' 

which has implicit dependence on T through Z{Xh), and is constant in T if and only 
if Z{X) is a constant in A, that is, if the axion does not couple to the dilaton. Indeed, 
as we mentioned at the end of section |^, the behavior of Z{X) as a function of T 
determines the behavior of Tcs/{sT /N^) as a function of T. In particular, on the 
large black hole branch, Tcs/isT/N^) is bounded from below by its value in the 
T — )■ oo limit, 

lim = (3.25) 

T^oo sT/N^ 27r ^ ^ 

If we use the preferred value k^Zq = 33.25 [42] then k,'^Zq/{2tt) ~ 5.29. Moreover, 
^cs/{sT/N^) will increase monotonically as T approahces Tc from above. 

For the simplest choice of Z{X), given in eq. (|2.12|) , Tcs/ {sT /N'^) has extremely 
mild dependence on T: as T approaches from above, Tcs/isT/N^) is nearly 
constant, experiencing an increase of only about 0.01%, mostly between 2Tc and T^, 
as shown in fig. ^. In bulk terms, the reason for this mild T dependence is that 
between T — t- oo and T = Tc, Xh increases from zero up to only Xh ~ 0.14, which for 
the Z{X) in eq. (p^.l2D translates into a very small change in Tcs/isT/N^). 



rcs/{K-Zo/2n) 
1.00020 h 



1.00015 - 



T/r, 



Figure 4: Our numerical result for Tcs / {sT /N^)-, normalized to the T — )• oo value 
k'^Zo/{2'k), as a function of T/T^ for the Z{\) given in eq. (|]l|), with a = 0.26 [42]. As 
T decreases, Tcs/isT/N^) remains nearly constant, experiencing only an approximately 
0.01% increase, mostly between 2Tc and Tc. 



On the other hand, for the Z{X) in eq. (|2.19| ), for different values of the coef- 



ficients Ci and C4 we find more variation in Tcs/isT/N^) as T approaches from 
above, as shown in fig. ^ For example, if Ci = 40 and C4 = 600, then Tcs/{sT/N^) 
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increases near Tc by more than a factor of six. For all values of ci and C4 that we con- 
sidered, most of the increase occurs between 2Tc and Tc. Fig. ^ shows Tcs/{sT/N^), 
normalized to the T ^ 00 value i^Z^j (27r), as a function of T/Tc for values of C\ and 
C4 that reproduce the lattice results for axial glueball mass ratios to within one sigma, 
eq. ( p.20| ). At the upper limits of the allowed (01,04) values, namely (ci, C4) = (5, 50), 
we find that as T approaches T^. from above, VcsK^T /N^) increases by about 60%, 
with most of the increase occuring between 2Tc and T^. 
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Figure 5: (a) Our numerical results for Tcs/isT/N^), normalized to the T — t- 00 value 
AC^Zo/27r, as functions of T/Tc, for the Z{X) in eq. ( |2.19D , for different choices of the 
dimensionless parameters (01,04). From the bottom (red) curve to the top (blue) curve, 
(ci,C4) = (0,0.26), (0.5, 0.87), (1,2.2), (5, 24), (10, 75), (20, 230), (40, 600). (b) Close-up of 
the curves for (from bottom to top) (01,04) = (0,0.26), (0.5,0.87), (1,2.2). In all of these 
cases, as T approaches from above Tcs/isT/N^) increases by anywhere from 0.01% up 
to a factor greater than six. The increase occurs mostly between 2Tc and Tg. 

In heavy ion collisions at RHIC and LHC, T reaches two to four times the QCD 
crossover temperature. We would thus like to know the value of Fes in QCD near 
the crossover temperature, which is a key ingredient determining the magnitude of 
any current produced via the CME [18].^ No controlled calculation of Fes from QCD 
at these temperatures exists, hence we turn to holography. Suppose we use A/" = 4 
SYM as a holographic proxy for QCD near the crossover temperature. The result 
for Fes in large- iVc, strongly-coupled M = 4 SYM is [30], 

^^s' = ^sT\ (3.26) 

Being a conformal field theory. A/" = 4 SYM has no phase transitions at non-zero T, 
so to obtain a sensible result we should consider the dimensionless quantity Tcs/T^- 
As a crude estimate we take as = / (47r) = 0.5 and we use Nc = 3, so that Xt = Gn, 
in which case we find 

F^g=7T^ ~ 0.045. {Xt = Gtt) (3.27) 



^We thank D. Kharzeev for a discussion on this point. 
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Figure 6: Our numerical results for Tcs/{sT/N^), normalized to the T 

^2 ' 



oo value 

K^Zo/27r, as functions of T/Tc, for the Z{X) in eq. ( ^.19 ) with (ci,C4) constrained such 



that the holographic model reproduces the lattice results for axial glueball mass ratios to 
within one sigma: < ci < 5 and 0.06 < C4 < 50. A generic choice of (ci, C4) within these 
limits will produce a curve inside the shaded region. The lower bound of the shaded region, 
given by the solid pink curve, has the lowest values, (01,04) = (0,0.06), while the upper 
bound, given by the solid blue curve, has the largest values, (01,04) = (5, 50). At the upper 
bound we see that as T approaches Tc from above, Tcs/{sT/N^) x (In) / {k^ Zq) increases 



by about 60%. The dashed line is the result for the optimal values (01,04) 
shown also in fig. ^ 



(0,0.26), as 



For a better estimate, let us consider Tcs(Tc)/T^ in IHQCD. As discussed above, if 
Z{X) is monotonic in A, then Tcs/isT/N^) is bounded from below by its value in 
the T — )• 00 limit, eq. ( p.25| ). We can obtain a lower bound on Tcs(Tc)/T^ by using 
the large-iVc YM lattice result for the entropy density at [47], s{Tc) = O.SIN^T^. 
Letting Ac denote the value of A^ at Tc, we find 



Tcs{Tc)/T^ = 0.31 X 



K^Z{X, 



> 0.31 X 



2tt 



1.64, 



(3.28) 



which is about 36 times larger than the A/" = 4 SYM estimate, eq. ( |3.27D . In fact, 
eq. ( |3.28| ) is closer to the perturbative QCD result, if we naively extrapolate to 
as = 0.5: rcs{T)/T^ ^ 30af ^ 0.94 (up to logarithms) [8-11]. If we consider the 
Z{X) in eq. (|2.19|) , and constrain ci and C4 to the values in eq. ( p.20|) , then we can 
also place an upper bound on Tcs(Tc)/T^, given by the solid blue curve in fig. For 
these choices of Z(A), we thus find 



1.64 < Tcs{T,)/T^ < 2.8. 



(3.29) 



Finally, we have also calculated Tcs using the small black hole solutions [39, 
40,42]. Our results for those cases appear in the appendix. Although the small 
black hole branch is always thermodynamically disfavored, we can actually use the 
results for rcs/('5T/A'^) on the small black hole branch to argue quite generally 
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that on the large black hole branch Tcs/isT/N^) should increase as T approaches 
Tc from above. A similar argument also applies for the bulk viscosity, as discussed 
in ref. [43]. The key result, shown in fig. |Tl] in the appendix, is that for T > 
Tcs/isT/N^) is larger on the small black hole branch than on the large black hole 
branch, but the two branches meet at Tmin- On the large black hole branch, then, 
Tcs/isT/N^) must increase as T — )■ T^m from above, in order to meet Tcs/isT/N^) 
from the small black hole branch. In fact, we can show in full generality that on 
the large black hole branch Tcs/{sT/N^) must increase as T — > Tmin- we simply 
take {d/dT)iTcs/isT/N^)) = {dXh/dT){d/dXh){K'^Z{\h) /2n) and observe that by 
definition {dXh/dT) diverges when T — )■ T^m, while {d/dXh){n'^Z{Xh)/2'7i) remains 
finite. Notice also that Tcs/isT/N^) itself remains finite when T — T^m. Given 
that Tmin is generally very close to we are then guaranteed that Tcs/{sT/N^) 
will be increasing as T Tc from above, if we assume that Z{X) is monotonic as 
a function of T between Tmin and Tc. In principle, Z{X) could exhibit maxima or 
minima for T G {Tmin,Tc), although such behavior seems un-natural. On the large 
black hole branch an increase of Tcs/isT/N^) as T Tc from above seems to be 
the generic behavior. We thus learn that the increase in Tcs/isT/N^) in the vicinity 
of Tc on the large black hole branch is tied to the existence of T^m, and hence to 
the existence of small black hole solutions. As argued in ref. [40], the existence of 
small black hole solutions follows from the fact that the zero-temperature theory is 
confining. These arguments suggest that perhaps any confining, strongly-interacting, 
large- Ai'c gauge theory with a (4+l)-dimensional holographic dual^*^ may exhibit an 
increase in Tcs/{sT/N^) in the vicinity of Tc. 



4. The Spectral Function 

We now turn our attention to Gr{u, k) with non-zero u and k. Generically Gr{uj, k) 
is a complex-valued function of the real variables u and k. A pole in Gji{u, k) 
indicates a large response to an infinitesimal source for q{x^), and is thus associated 
with a resonant excitation of the system. Being complex-valued, Gr{co, k) is not 
directly observable. To study the excitations of our system, we thus turn to the 
spectral function, —2\m.GR{(jj, k), which is real and hence observable in principle. 
Typically, a pole in Gr{uj, k) produces a peak in the spectral function. In this section 
we initiate the study of these peaks in our system. 

"'^'^Our arguments may not apply for (3+l)-dimensional confining theories obtained from higher- 
dimensional theories with compact spatial directions, such as the low-energy worldvolume theory 
on D4-branes with one spatial direction compactified and anti-periodic boundary conditions for 
fermions [29]. 

^^Given lr[\Gii{uj,k) we can obtain ReG/j(aj,fc) via a Kramers-Kroning relation, provided the 
large-w and large- |fc| asymptotics have been suitably regulated. 
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To be precise, we will compute Im Gr{u, k). To do so, we will compute Gr^u, k) 
using the membrane paradigm [32], as explained in section ^. In particular, we must 
solve eq. ( |3.21| ), which we reproduce here for convenience 



Z{\{r))h{rff\r 
with the boundary condition in eq. ( p.22D 



e + Z{\{T)fb{Tf\l-f{T)^' 



2 



(4.1) 



Cirn) = +lZ{\h)h{r^,)\ (4.2) 
and then obtain Gr{uj^ k) via eq. ( p.20|) . 



Gr{uj, k) = -K^Ml CO lim C(r, u, k). (4.3) 



We have not been able to solve eq. ([4.1|) exactly for all values of u and fc, hence we 



turn to numerical solutions. In this section we exclusively use the Z[X) in eq. ( |2.12| ), 
with C4 = 0.26. 

We consider first the case = 0. Fig. |^ shows our numerical result for Im Gr{u, k ■ 
0)/(TcMp) at Tc as a function of u/Tc. As we saw in section |^, for co sufficiently 
small, lmGR{u, k = 0) (x to. On the other hand, at asymptotically large to we expect 
lmGFi{u},k = 0) oc because in the UV the theory is conformally invariant and 
g(x'^) is dimension four. Our results are consistent with that expectation: fig. ^ shows 
that the function (1.6 x 10^^) x (w/Tc)^ *^^^ provides an excellent fit to our data. 

The scaling of lmGji{u, k = 0), and hence of Gr{co, k = 0), at asymptotically 
large a; is a divergence in the coincidence limit of the two-point function that prevents 
the correlator from obeying the sum rules and dispersion relations typically used to 
give physical meaning to the poles of Gr{u, k) in the complex cj plane, which require 

— * 

Gji{u!, k) to vanish at large frequency. Such a divergence may overwhelm peaks in 
lmGR{u, k), rendering them practically invisible. 

One way to improve the large-u; behavior of Im Gr{u, k) is to consider subtracted 
correlators. For example, one possible option is to determine the form of IthGrIoj, k) 
at large u exactly by solving eq. ([4.1|) in a WKB approximation, and then subtracting 
that large-w form from all subsequent calculations of Im G r{uj , k) . That approach 
encounters ambiguities in sub- leading divergences in u, as discussed for example in 
ref. [54]. We will instead eliminate the large-w divergence by computing Gr{u), k) at 
two different temperatures, Ti and T2, and then taking the difference, 

^Gr{uj, k- Ti, T2) = Gr{uj, k) - Gr{uj, k) . (4.4) 



T2 



Ti 



We could also imagine subtracting the T = result for Gr{ijJ, k), that is, by taking 
Ti = 0, but that is difficult to do numerically. When T = 0, Gr{u, k) is a sequence 
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Im Gr 




Figure 7: Our numerical results for lmGfi{uj, k = 0) / (TcMp) as a function oiu/Tc, at Tc, 
for the Z{\) in eq. with K^Zo = 33.25 and C4 = 0.26. The red dots are our numerical 

results while the solid blue curve is the function (1.6 x 10^'') x (u/Tc)^'^^^ . Our results are 
clearly consistent with the expectation that ImG/j(a;, A; = 0) oc o;"^ at large uj. 



of delta-functions whose locations and amplitudes correspond to the masses and 
wave-function normalizations of axial glueballs. We would need to subtract the 
enveloping function of this sequence of delta-functions, which is difficult to implement 
numerically. We will thus always consider Ti,T2 > Tc. Fig ^ shows our numerical 
results for lmGR{u,k = 0) at two different temperatures, Tc and 2Tc, while fig. || 
shows our numerical results for AlmGn^u, k = 0; Tc, 2Tc). In each figure we observe 
that the difference in ImGR^u, k = 0) between Tc and 2Tc approaches zero as co/Tc — )■ 
00, at least within our numerical precision. Our numerical subtraction thus appears 
to be reliable, so we may interpret peaks in lmGR{u, k) as physical excitations. 

From figs. ^ and we see that as T increases from Tc to 2Tc, \t[iGr{uj, A; = 0) 
changes by at most 10%. Fig. ^ also clearly reveals a minimum in AImG/j(co', k = 
0]Tc,2Tc) near u/Tc ~ 10 and a maximum near u/Tc ~ 22, indicating a shift in 
spectral weight towards higher w as T increases. Indeed, fig. ^ strongly suggests that 
a peak in the spectral function is moving to higher a; as T increases. The location 
of the peak, at u on the order of twenty times Tc, is roughly the same as the scale 
of the lightest O""*" glueball mass at T = 0, around 2600 MeV [53]. In other words, 
fig. ^ provides evidence that the plasma supports an excitation with roughly the same 
energy as the lightest O""*" glueball at T = 0. The width of the peak in fig. |]is about 
lOTc ~ 1300 MeV, so the excitation is reasonably long-lived. 



Figure 10 shows our result for the subtracted correlator with non-zero u and \k\ 



-22- 



Im Gr 



0.015 



T=2 r. 



0.010 



0.005 



>nii:::.........«»«"t* 



10 



15 



it) 



Figure 8: Our numerical results for ImG/j(a;, k = 0)/{TcMp) as a function of uj/T^, at Tc 
(lower blue dots) and at 2Tc (upper red dots), for the Z{X) in eq. (|2.12|) with k'^Zq = 33.25 
and C4 = 0.26. A both Tc and 2Tc, for lo/Tc sufficiently large ImG/j(a;, A; = 0) cx o;^. 

using the same two temperatures as above. We observe that as increases up to 
\k\/Tc ~ 10, the largest peak shifts from u/Tc ~ 22 up to u/Tc ~ 30. Although 
this change in the position of the peak is roughly order one, the change in the shape 
of the peak is very mild. In particular, the width of the peak changes very little, 
indicating that the lifetime of the excitation stays nearly constant as increases. 

The typical time scale for dynamical processes in the QGP created in heavy ion 
collisions is about Ifm/c ~ (200MeV)~^. Our results suggest the existence of a 
relatively long-lived excitation with energy on the order of 2600 MeV, corresponding 
to a time scale of about O.lfm/c. We cannot resist speculating that perhaps such 
an excitation, if present in the QGP, could dominate correlators of g(x^) and hence 
many dynamical CP-odd phenomena. Regrettably, we will leave a detailed analysis 
of this excitation, and its effect on CP-odd physics, for the future. 



5. Discussion and Outlook 

IHQCD is a state-of-the-art bottom-up holographic model for the low-energy physics 
of (3+l)-dimensional large- A^^c YM theory. In this paper we computed the retarded 
Green's function of the instanton density operator q{x'^) in the high-temperature, 
deconfined phase of IHQCD. Our primary motivation was to compute the Chern- 



- 23 - 



A(Im Gr) 




0;T„2T,)/{T,Ml) 



as 



Figure 9: Our numerical results for the difference AlmGji{uj,k 
a function of oo/Tc, for the Z(\) in eq. (|2l^ ) with k^Zq = 33.25 and C4 = 0.26. The 
difference goes to zero (within our numerical precision) as co/Tc — )■ 00, as expected. The 
prominent minimum at uj/Tc ~ 10 and maximum at w/Tc w 22 indicate a shift in spectral 
weight with increasing T, presumably from the motion of a peak in the spectral function. 




t igure 10: Our numerical results for AImGj?(w,A; = 0;Tc,2Tc)/{TcM^) as a function 
of uj/Tc and \k\/Tc, for the Z(A) in eq. with k^Zq = 33.25 and C4 = 0.26. As \k\ 

increases up to \k\/Tc ^ 10, the largest peak shifts from to /Tc ^ 22 up to w/Tc 30. The 
width of the peak changes very little. 



Simons diffusion rate, Fes, with the result in eq. ( |1.5|) . In particular, our result for 
Fes is proportional to Z{Xh), where Z{X) is the normalization factor of the bulk ax- 
ion action, and Xh is the value of the holographic 't Hooft coupling at the black hole 
horizon. A combination of available data for the topological susceptibility and axial 
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glueball spectrum of laige-Nc YM, and glueball universality, are sufficient to deter- 
mine the small and large A limits of Z{X) [37,38,41,44]. We considered several forms 
for Z{X). Assuming that Z{X) is a monotonic function of A, we found quite gener- 
ally that Tcs/{sT/N^) is bounded from below by its value in the T — )• oo limit and 
increases monotonically as T Tc from above. Indeed, we presented an argument 
that the same will be true in many (3+l)-dimensional, confining, strongly-coupled, 
large- iVc theories with holographic duals. For the Z{X) producing our optimal fit to 
the lattice results for the axial glueball spectrum, we found that the increase was 
only 0.01%. Fixing Z{X) completely by a least-squares fit to lattice results for the 
Euclidean two-point function of q{x^), as explained in section H, is an important task 
for the future. We also presented evidence for a relatively long-lived excitation in 
the system with energy roughly on the order of the mass at T = of the lightest 
0""*" glueball, which prompted us to speculate that perhaps such an excitation could 
dominate CP-odd phenomena in the QGP created in heavy ion collisions. 

IHQCD is dual to pure large-A^^^ YM, so an important goal for the future is 
to include the effects of quarks in the holographic calculation of Fes- Some key 
questions are how the quark mass and chiral symmetry breaking affect Fes- The 
axial and vector flavor U{1) currents are dual to two U{1) Maxwell fields in the bulk, 
and the quark mass operator is dual to a complex scalar field, a tachyon, that is 
bi-fundamental under these two gauge fields. In the bulk, the axion couples to the 
axial U{1) gauge field and the to phase of the tachyon, as explained in refs. [55,56]. A 
solution for the tachyon describing either a non-zero quark mass or chiral symmetry 
breaking can thus influence the axion and affect Fes- 

Introducing flavors fields would also enable us to compute holographically the 
current produced via the CME. A preliminary requirement is a bulk solution describ- 
ing a magnetic field and a net chirality. 

We plan to study these and other related issues in the future. 
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Appendix: The Small Black Hole Branch 



As discussed in section ^, when T > IHQCD admits two branches of black hole 
solutions, large black holes and small black holes [40] . In this appendix we compute 
Fes using the small black hole solutions. 

We can determine the dependence of Tqs on T in the large-T limit of the small 
black hole solutions as follows. For generality, we will consider a dilaton potential 
V{X) whose large-A asymptotic form is V{X) oc A^/^ (logA)'^, with P a non-negative 
real number. In the body of the paper we used P = 1/2. From fig. |l|, we observe 
that for the small black hole solutions, when T is large, is also large. When is 
large, Xh is also large, in which case we can approximate Z{Xh) ~ ZqC/^X^ and hence 
Tcs/{sT/Nl) ^ n^ZoCiXi/{2Ti). As shown in refs. [37,38,41,44], in the r ^ cx) 
limit, A(r) oc exp(r^/'^^~'^)) r * Evaluating at ry, gives us Xh in terms of r^^. From 
ref. [40] we know in terms of T on the small black hole branch in the r/^ — t- oo 
limit, Th oc T^^^^^^^ . We thus find 



cs 



sT/N- 



oc 



271 



\?,C(T/T,)T 



where C is a dimensionless positive constant that depends on the choice of V^(A). 



sT/N/ 



(a) 



(b) 



Figure 11: (a.) Our numerical result for Tcs/{sT/N^), divided by Zqk^ /{2tt), as a 
function of T/Tc, for the Z{X) in eq. with C4 = 0.26. The upper dot-dashed blue 

curve is our result obtained from small black hole solutions, while the lower solid blue 
curve is the result obtained from large black hole solutions. Both curves begin at Tmin, 
indicated by the vertical dashed black line, which is slightly below Tc. The result on the 
small black hole branch increases as T increases, and in the T — )• oo limit approaches the 
form in eq. (|l]). (b.) Close-up of (a.) near Tmin- 
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To compute Fes in the entire range Tmin < T < cxo, we resorted to numerics. 
For the Z{X) in eq. (|2.12|) with C4 = 0.26, our results appear in fig. |Tl|, where we see 
clearly that the result grows as T increases, and in the T — )■ cxd limit approaches the 
form in eq. Fig. 11 also shows that the result for Tcs/i^T/N^) is always larger 
on the small black hole branch than on the large black hole branch, except at T^m 
where the two are equal. This result is important for our argument at the end of 
section ^ that Tcs/isT/N^) computed on the large black hole branch will increase 
as T approaches Tc from above. 
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